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Q . ABSTRACT 

' Using a statistically representative sample of 911 central galaxies (CENs) from the SDSS 

• , DR4 Group Catalogue, we study how the structure (shape and size) of the first rank (by stellar 

• mass) group and cluster members depends on (1) galaxy stellar mass (Mgtar), (2) the global 

environment defined by the dark matter halo mass (Mhaio) of the host group, and (3) the lo- 
Q . cal environment defined by their special halo-centric position. We establish a GALFIT-based 

5— ( ■ pipeline for 2D Sersic fitting of SDSS data to measure the Sersic index, n, and half-light ra- 

C/3 I dius, r^o, from r-band galaxy images. Through tests with simulated and real image data, we 

. demonstrate that our pipeline can recover galaxy properties without significant bias. We also 

find that uncertainties in the background sky level translate into a strong covariance between 
[ the total magnitude, the half-light radius, and the Sersic index, especially for bright/massive 

^ ■ galaxies. We apply our pipeline to the CEN sample and find that the Sersic index n of CENs 

' depends strongly on Mgtar, but only weakly or not at all on Mhaio- The n-Mgtar relation holds 

Vf^ , for CENs over the full range of halo masses that we consider Less massive CENs tend to be 

disk-like and high-mass systems are typically spheroids, with a considerable scatter in n at 
all galaxy masses. Similarly, CEN sizes depend on galaxy stellar mass and luminosity, with 
early and late-type galaxies exhibiting different slopes for the size-luminosity (t^q-L) and the 
' size-stellar mass (rsp-Mstar) scaling relations. Moreover, to test the impact of local environ- 

Q^ I ment on CENs, we compare the structure of CENs with that of comparable satellite galaxies 

■ (SAT). We find that low mass (< lO^^'^^h^^M©) SATs have somewhat larger median Sersic 

indices compared with CENs of a similar stellar mass. Also, low mass, late-type SATs are 
moderately smaller in size than late-type CENs of the same stellar mass. However, we find 
no size differences between early-type CENs and SATs and no structural differences between 
i-j I CENs and SATs when they are matched in both optical colour and stellar mass. The similarity 

in the structure of massive SATs and CENs demonstrates that this distinction has no signif- 
icant impact on the structure of spheroids. We conclude that Mgtar is the most fundamental 
property determining the basic structural shape and size of a galaxy. In contrast, the lack of a 
significant rt-Mhaio relation rules out a clear distinct group mass for producing spheroids, and 
the morphological transformation processes that produce spheroids must occur at the centres 
of groups spanning a wide range of masses. 

Key words: galaxies: clusters: general — galaxies: evolution — galaxies: formation — galax- 
ies: fundamental parameters — galaxies: structure. 



1 INTRODUCTION 
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Understanding the role that environment plays in shaping the mor- 
phology of galaxies remains an important challenge in the study of 
galaxy formation and evolution. In the standard galaxy formation 
and evolution paradigm, all galaxies started as star-forming disks 
at the centres of small dark matter (DM) haloes. Subsequent hierar- 
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chical evolution has transformed galaxies into spheroids to produce 
the bimodality observed in the present-day populations. Tied to 
the strong colour bimodality are differences in morphology and re- 
lated galactic structure. Blue star-forming galaxies tend to be disk- 
dominated with exponential radial light profiles (late-types), while 
red non-star-forming systems have spheroidal mass distributions 
with steeper light profiles (early- types). Many observations suggest 
that there may be an environmental component to differences in 
galaxy properties; e.g., dependence of morphology (e.g.. Dresslen 
Il98d iGoto et alj|2003l ; iMcIntosh et alJl20Cl4 iBlanton et al.| 2005 ) 
and star formation related properties (e.g., iHashimoto et al.lll998 



iBalogh et all |2004 . iKauffmann et al.l |2004|) to local galaxv den- 
sity. Presumably, the physical processes that have built larger DM 
haloes housing groups and clusters of galaxies must be responsible 
at some level for the observed evolution of the two primary galaxy 
populations. 

The latest theoretical models of galaxy evolution invoke the 
transformation of blue disks into red spheroids within a hierar- 
chical framework to explain the factor of two growth observed 
in the early-type galaxy (ETG) population since z = 1_ (e.g 



ixy , , 

Bell et al.'20Q4VBla ntonl200^ : lBorch et alj2006l ; lFaber et al.l200'7l : 

Brown et al. 2007). It is usually assumed that to successfully re- 
produce the galaxy bimodality requires physical mechanisms that 
both transform star formation and morphology. A variety of galaxy 
transformation scenarios can be found in the literature, many of 
which are predicted to be important only in particular environ- 
ments. Yet, recent advances in understanding bimodality have fo- 
cused on environ mental processes that mostl y impact star forma- 
tion. For example. lvan den Bosch et al. I( l2008l) successfully demon- 
strate that quenching is important for producing redder galaxies 
over a range of different environments, but a clear picture of what 
governs morphological bimodality is still lacking and controversial. 
To shed light on whether a special environment exists for the trans- 
formation of galaxy morphology, we study the shapes and sizes of 
a representative sample of central galaxies from galaxy groups an d 
clusters in the Sloan Digital Sky Survey (SDSS. lYork et al.ll2000h . 
Quantifying the structural properties of galaxies, such as the steep- 
ness of the light profile shape, provides a direct means to assess 
morphological transformation. 

If all galaxies started as small disks, then the simple existence 
of spheroidal systems makes it clear that morphological transfor- 
mation occurs. The transformation from late to early-type galaxies 
may occur in one traumatic event or be the result of multiple pro- 
cesses over billions of years. There is no shortage of theoretical 
predictions for the creation of ETGs, from violent galaxy-galaxy 
mergers to the slow fading of disks. Generally speaking, the early 
type population includes a range of spheroid-dominated morpholo- 
gies with and without a disk component (e.g., ellipticals, lenticu- 
lars, and bulge-dominated spiral galaxies). Thus, it is likely that 
more than one physical mechanism is responsible for the variety 
among ETGs. Narrowing down the processes that are most respon- 
sible for turning disks into spheroids remains a critical piece of the 
galaxy evolution puzzle. In this paper, we approach this problem 
by testing whether or not there is a specific environment where a 
strong morphological transition takes place that can then be tied to 
a particular process. 

There are several morphologically-altering mechanisms that 
are predicted to be effective mainly in the high-density environs 
of massive groups and clusters. Harassment, the cum ulative effect 
of ma ny high speed encounters with satellite galaxies jMoore et al.l 
is predicted to occur primarily in groups and clusters after 
a satellite galaxy is accreted, and may transform a disk galaxy into 



a more early-type morphology by heating and 'puffing up' the disk 
component. Tidal stripping, the effect of the tidal force suffered 
by a satellite along its orbit, may transform a satellite galaxy into a 
spheroid by removing its disk, and it may be effective in haloes over 
a large mass range. Since the above two processes change the stellar 
mass of a galaxy by at most a factor of two, if they are responsible 
for the morphological transformation one would expect to see satel- 
lites with statistically significant differences in their morphological 
and structural pr operties compared wi t h cent rals o f a similar stellar 
mass. However, Ivan den Bosch et al] ( 1200 8l) and IWeinmann et al.l 
1 I2OO8I) find that satellite galaxies are only marginally more concen- 
trated than central galaxies with same stellar mass, which suggests 
that satellite galaxies only undergo a minor change in their mor- 
phology after they fall into a massive cluster, and that the above two 
processes may not be sufficient to explaining their major morpho- 
logical transformation. It is important to note that these processes 
basically produce early types by diminishing the disk structure of 
later-type spirals; i.e., these mechanisms may produce SO/Sa early 
types, but they do not create massive spheroids or elliptical galax- 
ies. 

Besides the aforementioned cluster-specific mechanisms that 
appear to have only minor impact on the morphologies of satel- 
lite galaxies, numerical simulations demonstrate that the merger of 
two d isk galaxies wi t h similar masses produces a spheroid galaxy 



'2003!; 



Toomrd I1977I; [Barnes & Hernguistl Il996l ; iNaab & Burkert 



Cox et al 



200fT) and, hence, are likely an important mech- 
anism for the formation of spheroids and ellipticals. Major merg- 
ers are believed to be efficient in group-size haloes and to be sup- 
pressed in more massive haloes because of the increasing differ- 
ences in the relative velocities between member galaxies com- 
pared to their own internal velocity dispersions. It has long been 
assumed that the smaller velocity di spersions found in g alaxy 
groups allow more galaxy interactions jCavaliere et alj[l992l) ; also 
the orbital decay timescale is shorter in l ower-mass halo e s (e.g . 
ICoorav & Milosavlievid |2005|) . However, iMcIntosh et"aLl ( 12003) 
find that major mergers among massive galaxies occurs at the cen- 
tres of clusters, as well as large groups, yet this merely makes 
bigger spheroids from smaller spheroids. Moreover. iHopkins et al.l 
also show, using theoretical arguments, that the massive 
central galaxies of clusters still have a large chance to merge with 
their satellites. These findings highlight once again the open ques- 
tion of which environment or halo mass is ideal for transforming 
disks into spheroids. It appears that major mergers preferentially 
happen between central and satellite galaxies; therefore, the exis- 
tence of a special halo mass for major morphological transforma- 
tions such as mergers might manifest itself as a noticeable change 
in the morphological distribution of central galaxies at some spe- 
cific halo mass. We already know that central galaxies living in 
small haloes have disk-like shapes with low concentrations and flat 
light profiles while those in large haloes have spheroid-like shapes 
with high concentrations and steep light profiles. A careful study of 
the distribution of structural properties in haloes spanning a range 
in mass will help shed light on this open question. 

Our analysis is based on two quantitative measures of galac- 
tic structure, the Sersic index and the size, which are directly 
related to galaxy morphology. The measureme nt of galaxy light 
profile shapes and sizes has a long hi story ( Shaw & Gilmord 
1989 ; Byun & Freeman 1995; Ide Jond 1996; Simard 199|; 



Khosroshahi et al.[|20od and references therein) and the develop- 



ment of automatic routines to handle the huge number of galax- 
ies from modern surveys is well-motivated. In this work, we de- 
velop a powerful pipeline for applying a well-tested and popular 
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software package for two-d imensional (2D) galaxy image fitting 
(GAL FIT.IPeng et al.ll2002l) to SDSS r-band data. We fit a Sersic 
model l lSersidl968l) tothe images of galaxies and use the Sersic in- 
dex, n, and half-light radius, 7-50, of the best fit models to describe 
the shapes and sizes of galaxies, respectively. A well-known, yet 
often overlooked, issue in profile fitting is the critical sensitivity t o 
the estimate of the background sky level jMacArthur et alj[2003h . 
We explore this issue in detail and discuss its impact on quantify- 
ing the structure of galaxies, in particular that of high-mass central 
galaxies. 

In this paper, we try to answer two questions about the en- 
vironment of morphological transformation: (1) is there a critical 
DM halo mass where central galaxies are transformed from late- 
type to early-type? (2) is the central position in groups and clusters 
a special place for determining the structure of galaxies. Previ ous 
work arguing wheth er galaxy morphol ogy (e.g., Dressledl 19801) or 
star formation (e.g. .I van der Welllioosh depends more critically on 
environment employ local galaxy density measurements that are 
less natural and less physically meaningful than the halo mass and 
the location within a host halo jWeinmann et al. 2006t). H ere we 
use the SDSS DR4 Group Catalogue of l Yang et all ilOOHt . which 
provides statistical measures of the host halo mass (global environ- 
ment) and the halo-centric position (local environment) for SDSS 
galaxies. We study the quantitative structure of galaxies living at the 
presumed dynamical centres of DM haloes spanning nearly three 
orders of magnitude in mass. In hierarchical models of galaxy for- 
mation, all galaxies begin as the central galaxy in a smaller halo and 
then become satellite galaxies if the small halo merges with a big- 
ger one. Therefore, comparing the structural properties of centrals 
and satellites can shed light on the evolution of galactic morphol- 
ogy and its dependence on environment, and help answer whether 
central galaxies are a distinct population with unique formation his- 
tories. 

In !j2]we present our sample selection. In ^we present and 
test our fitting pipeline and show the sensitivity of the fits to the 
background sky level. We present the results of our fits to the cen- 
tral galaxies in ^ both their Sersic indices and sizes, and compare 
them to satellite galaxies. We summarise our main conclusions in 
^ We also include an appendix where we compare our fits to those 
of the NYU-VAGC (Blanton et al. 2005). Throughout we adopt a 
flat ACDM cosmology with Qrn = 0.3, f^A = 0.7 and use the 
Hubble constant in terms of h = _ffo /100km Mpc~^. 



2 SAMPLE SELECTION 

2.1 Central Galaxy Sample from SDSS Group Catalog 

To study the structural properties of central galaxies (CENs) 
in group s and clus t ers, w e use the SDSS galaxy group cata- 
logue of I Yang et alj ( l2007l hereafter Y07). This catalogue is con- 
structed using the Ne w York University Value-Added Galaxy 
Catalog (NYU-VAGC, iBlanton et a l7 '2005^ reprocessing of the 
spectroscopic 'Main' selection ([strauss et al. 2002) from the 
fourth data release of Sloan Digital Sky Survey (SDSS DR4, 
lAdelman-McCarthv et al.1l2006i) . The NYU-VAGC provides im- 
proved reductions and additional galaxy property measurements. 
The group catalogue provides two physically motivated measures 
of environment for each galaxy: (i) the dark matter halo mass of its 
host group (Mhaio), and (ii) its distance to the central (i.e, highest 
stellar mass) group member of the group. Here we briefly describe 
the group catalogue and we discuss its redshift completeness as it 
relates to our final selection of a representative sample of CENs. 



The details of the group finder used to construc t the DR4 ver - 
sion of the group catalogue are given in I Yang et al] ilOOj. \200% . 
Briefly, the group finder starts with a friends-of-friends (FOP) al- 
gorithm using a restrictively small linking length to define poten- 
tial groups and their centres in a galaxy redshift survey. A rough 
group mass is estimated from the total group luminosity, assuming 
a mass-to-light raticQ. From the group mass, the group-finder uses 
an adaptive filter to iterate on the virial search parameters (pro- 
jected radius and velocity dispersion), which in turn are used to se- 
lect group members in redshift space. This method is iterated until 
the group members converge. This algorithm has been thoroughly 
tested with mock galaxy redshift surveys and is shown to be more 
successful than conventional FOP finders. The average complete- 
ness of individual groups in terms of membership is ~ 90 percent, 
with only ~ 20 percent contamination from interlopers. The halo 
mass of each group is estimated using two methods: (1) the ranking 
of its total characteristic luminosity and (2) the ranking of its total 
characteristic stellar mass. As shown in Y07, both methods agree 
very well with each other, with a scatter of 0. 1 dex (0.05 dex) at low 
(high) halo masses. In this paper, we use the halo mass from stellar 
mass ranking. Finally, owing to the r < 17.7 (extinction-corrected) 
magnitude limit of the Main sample, the minimum Mhaio for which 
the group selection is complete changes with redshift. Therefore, 
we limit our analysis to groups with z < 0.08, which are detected 
with high completeness down to log(Mhaio/h~^M0) — 11.78, 
allowing us to study galaxies from small groups with good image 
resolution. 

The spectroscopic completeness of the sample used to deter- 
mine galaxy groups plays a crucial role in identifying CENs. Owing 
to fibre collisions, the Main sample misses about 8 pe rcent of SDSS 
galaxies meeting the spectroscopic targeting criteria iBlanton et al.l 
l2003i). This effect is se vere in regions of high galaxy number den- 
sity faogg et al.ll2004h , such as in large groups and clusters. The 
Y07 group catalogue contains three samples that address this in- 
completeness differently. Each sample spans the redshift range of 
0.01 < z < 0.20. Sample I contains 362356 Main galaxies. Sam- 
ple II includes 7091 additional galaxies wit h spectroscopic red - 
shifts measured from other surveys (e.g., 2dF. IColless et aklboOll) . 
Sample III adds 38672 galaxies missing redshifts due to fibre- 
collisions, which are assigned a redshift based on their nearest 
neighbour For the three samples, the group finder detects 201621, 
204813 and 205846 groups, respectively. 

Note that the fibre-collision correction applied to sample III 
may transfer some CENs in sample II to satellite galaxies in sam- 
ple III and vice versa. For our selection, we choose only galaxies 
identified as CENs in both sample II and sample III. These galaxies 
are guaranteed to be CENs regardless of whether the fibre-collision 
correction was applied or not. In addition, we use the Mhaio of 
these CENs drawn from sample II. Our analysis (discussed in Sj3]l 
is based on a CPU-intensive galaxy image fitting routine, thus, to 
construct a representative sample we randomly select CENs from 
halo mass bins spanning the full range of groups in the z < 0.08 
volume-limited sample. For log(Mhaio/h"^Mo) = [12.0 : 14.0], 
we randomly select 100 CENs from eight bins of 0.25 dex width. At 
higher halo masses the number of groups decreases rapidly, thus to 
maintain good statistics we select 100 CENs from the [14.0:14.5] 
bin, and use all 11 CENs in groups with log(Mhaio/h~^MQ) > 



^ The resuhing group catalogue is insensitive to the initial assumption re- 
garding the M/L ratios. 
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Figure 1. The colour-Mstar and Mstar-Mjjaio distributions of our selected CEN sample (black circles). Contours show the distribution of all z < 0.08 CENs 
from the SDSS DR4 group catalogue. Contours with different colours connect bins with same number of galaxies as labelled. Each bin has a width of 0.03 in 
colour and 0.05 dex in Mgtar in the left panel, and a width of 0.05 dex in Mgtar and 0.06 dex in MjjjjIq in the right panel. 



14.5. The total number of CENs in our representative sample is 
911. 

In Figure [T] we plot the stellar mass (Mstar) distributions of 
our CEN sample as a function of galaxy colour and halo mass. For 
each galaxy, the NYU-VAGC provides the Petrosian "'^(g — r) 
colour, K+E corrected to 2: = 0.1 (see Y07 for d etails), and we es - 
timate Mstar using the colour-derived M/L from lBelfet aLl( l2003h . 
As shown in Figure[T] our representative CEN selection samples the 
same Mstar, Mhaio and "'^ (<? — r) space as all CENs with z < 0.08 
in the group catalogue (shown by the contours). We note that our 
selection of CENs spans a wide range of Mstar running from 10^ * 
to lO^^ '^h'^MQ. Owing to the exponentially decreasing number 
density of high-mass haloes and the Mstar-Mhaio relation of CENs 
(Fig. 1, right panel), our selection of a constant number of groups 
per Mhaio bin results in a CEN sample that is heavily-weighted 
towards higher stellar masses and redder colours. 

2.2 Matched Satellite Galaxy Samples 

In addition to studying the structural properties of CENs, we 
also want to investigate whether their central position in groups 
or clusters produces a distinctive structural difference com- 
pared to non-CEN galaxie s. Following a similar method as in 
Ivan den Bosch et al] | |2008 !). we construct two control samples of 
satellite (SAT) galaxies, one to match CENs in stellar mass only, 
and the other to match in both stellar mass and colour. For each 
CEN in our sample, we first randomly select from all SATs with 
z < 0.08 a similar-mass counterpart with a Mstar within ±0.08 
dex, hereafter the SAT sample SI. We next find SATs likewise 
again matched to CENs in Mstar, but further restricted to match 
their colour within A"'^ — r) = ±0.03, hereafter the SAT sam- 
ple S2. Once a SAT is matched to a CEN, we remove it from the 
pool so that there is no duplication in each SAT sample. The match- 
ing criteria are equal to the measurement uncertainties of the stellar 



mass and the colour ( iBell et alj|2003l) , and using a matching crite- 
ria with either a smaller AMstar or dsP'^{g ~ r) will greatly reduce 
the number of matched SATs, especially towards the massive end, 
because we use a volume-limited sample and don't allow duplicate 
matches. For CENs with log(Mhaio/h"^MQ) < 11.0 we find a 
matched SAT using the above criteria more than 90 percent of the 
time. At higher masses (log(Mhaio/h-'MQ) > 11.3), SATs be- 
come rare and the fraction of massive CENs with a matched SAT 
rapidly drops to less than 10 percent. Hence, we achieve SAT sam- 
ples of 769 (SI) and 746 (S2) matches. 

In Figure |2l we plot the stellar mass and colour comparisons 
for our CEN-SAT matched samples as a function of CEN stellar 
mass. There is a slight bias in that SATs in both sample S 1 and S2 
are less massive than their counterpart CENs with a median dif- 
ference (SAT - CEN) of -0.01 dex, and SATs in SI are obviously 
redder than their counterpart CENs. The difference of colour for 
the SAT sample S2 {lower right panel) is also small, with an almost 
zero median for low mass (log(Mstar/h"^M0) < 10.5) SAT-CEN 
pairs and a median of 0.005 for massive (log(Mstar/h-2M0) > 
10.7) pairs. In the lower left panel we see an obvious difference be- 
tween the colour of SAT-CEN pairs in SAT sample SI, with a me- 
dian difference of 0.03 for massive (log(Mstar/h-2MQ > 10.7) 
pairs and 0.3 for low mass (log(Mstar/h~^MQ) < 10.3) pairs. 

To evaluate whether the above biases are statistically signif- 
icant or not, we match our CEN sample with mock samples cre- 
ated from itself with random changes on Mstar and colour accord- 
ing to the measurement uncertainties. The match is done with the 
same criteria used for constructing the SAT sample S2 and re- 
peated 20 times to produce a distribution of medians of AMstar 
and A"'^ — r) per Mstar bin between the CEN sample and mock 
samples. We plot the results in Figure|2]with circles and errorbars 
showing the mean and 3a deviation respectively. We can see that 
the biases of our real CEN-SAT matches are within 3cr deviations 
of the self-matching results, except in the lower left panel, where 
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Figure 2. The difference of Mstar {upper panels) and colour {lower panels) between matched CENs and SATs as a function of CEN Mstar- The left column 
shows the differences for the SAT sample SI, which is matched with our CEN sample only in stellar mass, while the right column shows the SAT sample S2, 
which is matched with our CEN sample in both stellar mass and colour. The solid line in each panel shows the median of the difference for each 0.25 dex wide 
Mstar bin. For comparison, circles with errorbars show the mean and 3cr deviation of the biases in each mass bin from matching the CEN sample to 20 mock 
samples created from itself with random changes according to measurement uncertainties (see text for details). 



the bias at log(Mstar/h~^MQ) < 11 are much larger than the 3a 
deviations. This result suggests that the Mstar and colour distribu- 
tions of our SAT samples are not significantly different from the 
distributions of our CEN sample, except for the case of colour in 
our SAT sample S 1 , where we find that SATs are redder than their 
CEN counterparts as Ivan den Bosch et"^ j2008l) finds. However, 
the colour difference between low mass CENs and SATs in sample 
SI (A.^'^{g — r) ~ 0.3) is larger than the average difference of 
0.1 in the same stellar mass region found bv Ivan den Bosch et al.l 
( l2008h . The reason for this disagreement could be the small num- 
ber statistics in our SAT sample. We conclude that the SATs from 
sample S2 is effectively matched in Mstar and colour with our CEN 



sample, while the SATs from sample S 1 are only matched in Mst^ 
and have a systematic redder colour than the CEN sample. 



3 MEASURING THE STRUCTURE OF GALAXIES 

The primary aim of this paper is to quantify the structural properties 
of CENs. In particular, we measure the shape an d size of the 2-D 
luminosity profile of each galaxy using GALFIT aPeng et alj2002h . 
This code fits a parametric model to the surface brightness profile of 
a galaxy image and o utputs a set of best-fitting parameters. For our 
analysis we adopt the lSersid ( Il968h model to describe the surface 
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brightness at radius r of a galaxy, E(r) = EeC"''"'''''''^'^''' '' 
where rso is the half-light radius of a galaxy, Ee is the surface 
brightness at rso, n is the Sersic index, and &„ is coupled to 
n such that half of the total flux of a galaxy is within rso (for 
0.5 < n < 10, bn ~ 2rz — 0.327). In addition to rso and n, 
GALFIT outputs the best-fitting total magnitude mtot, axis ratio 
h/a, and position angle 6. The Sersic profile is routinely used for 
galaxy structure analysis to provide the half-light radius and the 
Sersic index measures the galactic light profile shape; e.g., n = 4 
is the de Vaucouleurs' r^''* profile and n = 1 is the exponential 
disk profile. We choose GALFIT as our fitting tool because it can 
simultaneously fit Sersic profiles to several galaxies in one image, 
which is advantageous for galaxies in dense environments, where 
galaxies have a high chance to be overlapping with one other. In 
what follows, we outline our image fitting pipeline, test it with sim- 
ulated galaxies, discuss technical issues, and estimate the effects 
that background sky estimation has on parameter uncertainties. 



3.1 The Galaxy Fitting Pipeline: Modified GALAPAGOS 

To run GALFIT on each galaxy in our sample, we require a postage 
stamp image with an appropriate size to measure structure over the 
full extent of an object, the point spread function (PSF), the initial 
guess for the fitting parameters, and an estimate of the background 
sky level. In our pipeline, we start from the fully-processed SDSS 
imaging frames, which have a size of 2048 x 1042 pixels. We em- 
ploy GALAPAGOS (Barden in prep.) to process the whole image 
frame and to provide the needed information to GALFIT. GALA- 
PAGOS was originally designed to facilitate fitting large galaxy 
data sets based on HST/ACS images. We have modified this rou- 
tine to work on SDSS images. In brief, GALAPAGOS takes the 
following steps: (1) it detects neighbouring sources and produces 
image masks; (2) it cuts out postage stamps for detected sources; 
(3) it prepares an input file for GALFIT; and (4) it estimates local 
sky values for target galaxies. We outline the details of these steps 
below. 

Source Detection and Initial Fi t Parameter Guesses: We use 
SExtractor (Bertin & Arnoutj [T996h to detect and mask nearby 
companion sources in the SDSS image of each galaxy that we de- 
sire to fit, and to provide initial fit parameter guesses for the pri- 
mary galaxy and any close neighbours that will be simultaneously 
fit. SExtractor provides useful estimates of galaxy properties (mag- 
nitude, size, axis ratio, position angle) for calculating the initial 
guesses for the GALFIT parameters. A set of configuration param- 
eters defines how SExtractor detects sources. After tuning, we find 
that for SDSS r-band images the following configuration works 
best: DETECT JV[INAREA=25, DETECT.THRESH=3.0 and DE- 
BLENDJVIINCONT=0.003. This configuration provides a good 
trade off between detecting and deblending most bright and ex- 
tended sources without artificially deblending galaxies with strong 
substructures into multiple sources. Since our goal is to study bright 
galaxies in groups or clusters, we keep this configuration to al- 
low a high success rate on bright galaxies. To distinguish whether 
companions are stars or galaxies, all sources with SExtractor flag 
STAR_CLASS<0.9 are classified as galaxies. However, a small 
fraction of sources may be misclassified using this criterion, and 
we will discuss the effect of this misclassification on our fitting re- 
sults in ^331 

Postage Stamp Images: GALAPAGOS produces a rectangular 
postage stamp centred on each galaxy of interest. The purpose of 
postage stamps is to reduce the CPU time for fitting one galaxy. The 



postage stamps are cut to a size that is large enough to ensure that 
the outer light profile will be fit. GALAPAGOS uses parameters in 
the SExtractor catalogue to determine the X and Y dimension in 
pixels of the postage stamps for each object: 

X = 2.5 {\sin{e)\ + {l-e)\cos{e)\), (1) 

Y = 2.5 X a X rKron X {\cos{e)\ + (1 - e)\sin{e)\), (2) 

where a x rxron is the Rron radius along the major axis in units 
of pixels, 9 is the position angle, and e is the ellipticity. 

GALFIT Inputs: For running GALFIT, GALAPAGOS pro- 
duces an input file of initial parameter guesses for the fit- 
ting parameters based on the SExtractor output as follows: 
nitot.i = MAG_BEST for the apparent r-band magnitude, rso,i = 
0.162r}if^, where r/i„^ is FLUXJIADIUS; b/a- = 1 - e, where 
e is ELLIPTICITY; and Oi = THETA.IMAGE. We start with 
an initial Sersic index of n = 1.5. We note that for this analy- 
sis, we do not use the higher fitting modes, such as diskyness and 
boxiness, offered with the GALFIT software. Nearby companions 
within 1.5 times the SExtractor Kron aperture of the target galaxy 
are fit simultaneously with a Sersic model using initial parameters 
also determined as described above. Companions further away are 
masked out using the masks provided by SExtractor. In addition to 
the input file, GALFIT requires a PSF image to convolve with the 
model image. The PSF at the centre of the target galaxy is extracted 
from the SDSS photo pipeline by employing a SDSS published tool 
readAtlasImagefl 

Background Sky Estimation: The background sky level is a 
critically important ingredient in galaxy image fitting. For example, 
an overestimation of the sky can result in flux, size, and Sersic in- 
dex underestimation. GALAPAGOS includes a sophisticated way 
to measure the local sky around a g alaxy, which is demo nstrated 
to be successful for ACS images t Haussler et al. 1 12007|) . How- 
ever, GALAPAGOS uses a hierarchical iteration for fitting galax- 
ies, from bright to faint, over the whole frame to isolate sky pix- 
els. Rather than using this CPU-intensive approach, we rely on 
the SDSS, which provides useful and well-tested sky estimates. 
The SDSS global sky is considered to be a good measure of the 
background sky level for st udying the structure of SDSS galaxies 
Jvon der Linden et alj2007l) . We treat the background sky as a fixed 
flux pedestal during the fitting to reduce degeneracies between the 
sky and the outer isophotes of high-n models. 

3.2 Testing the PipeUne with Simulations 

We test our SDSS image fitting pipeline by running it on 850 sim- 
ulated galaxies. The goals of this test are two fold: (1) to estimate 
random and bias of the structural parameters returned by our fit- 
ting pipeline; and (2) to confirm that using the SDSS global sky 
does not produce bias. We compare the actual properties that de- 
fine each simulated galaxy (i nput) to the correspo nding fit result 
(output) following Figure 9 in lBlanton et al.l ( |2005|) . 

For our tests, we use the SDSS-based simulations of 
iBlanton et al] ( l2005h . Briefly, each simulated galaxy is an axisym- 
metric Sersic model. The simulation sample has a range of input 
parameters matching the NYU-VAGC Large-Scale Structure (LSS) 
sample from SDSS DR4. Each simulated galaxy is converted to raw 
data units, scaled to the SDSS pixel size, and convolved with the 
PSF at its position. After adding Poisson noise, the resulting image 

^ http://www.sdss.org/dr4/products/images/read_psf.html 
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Figure 3. A comparison between input and output parameters when using GALHT to fit Sersic models to 850 simulated Sersic galaxies. We plot the difference 
(output-input) as a function of input for three important Sersic parameters: the index n, the half-light radius rso, and the total magnitude. In each panel, grey 
points represent individual simulated galaxies. The solid lines show the 3rd quartile, median and 1st quartile (from top to bottom). 



is placed into an actual SDSS image at a random location. Gener- 
ally, about 60 simulated galaxies were added to each 2048 x 1024 
pixel image. 

The results of applying our pipeline to the SDSS simulations 
are shown in Figure[3] We plot the output-input offset as a function 
of input for three important parameters: the total magnitude mtot, 
the half-light radius rso, and the Sersic index n. The results demon- 
strate that our pipeline successfully recovers these structural param- 
eters of simulated Sersic galaxies with almost no offset and only a 



small scatter. For half of the simulation sample the pipeline returns 
structural measurements within 10 percent of the true value. As ex- 
pected, there is a larger parameter offset scatter for simulat ed galax- 
ies wi th Uin > 4 than for those with ni„ < 4 (see Haussler et al.l 
ioOTi)). We see a slightly increased scatter with magnitude owing 
to the lower S/N of fainter galaxies. We note that the good agree- 
ment between input and output parameters also demonstrates that 
the SDSS global sky is a good choice for galaxy image fitting. We 
address the effects of sky uncertainties on our structural properties 
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measurements in more detail in i]3.4| To the degree that Sersic pro- 
files provide a reasonable fit to the true light distribution of galax- 
ies, the good performance of our pipeline on simulated galaxies is 
promising for the analysis of real galaxies. 



3.3 Sersic Fitting of the CEN Sample 

We apply our pipeline to fit each galaxy in the CEN sample with 
a Sersic model. To evaluate the fit quality we visually checked a 
random subsample of 200 CENs and find that about 15 percent of 
our fits suffer from a variety of technical issues. These issues in- 
clude stars misclassified as galaxies and unreliable fits to compan- 
ion sources. Here we discuss what portion of our fitting results are 
affected by these issues and how we can correct for them. 

The accuracy of star/galaxy classification affects the fitting 
quality. Although CEN galaxies can always be correctly classified, 
a small portion of their companions may suffer from misclassifi- 
cation. Some galaxies are misclassified as stars and some stars are 
misclassified as galaxies, but only the latter has a severe impact 
on our fitting. Fitting a Sersic model to a stellar profile typically 
results in a very large n with an overestimated extent and an over- 
estimated flux. As a result, this false extension takes light away 
from the whole image and hence causes an underestimate of mtot , 
rso and n for the target galaxy. Conversely, companion galaxies 
misidentified as stars and fit by a PSF are spatially compact and 
faint, thus their improper treatment has little impact on the mea- 
sured structure or flux of the target galaxy. We estimate that less 
than 5 percent of our CEN sample (in a random sense) suffer from 
the issue of companion stars being misclassified as galaxies. 

A drawback to fitting several sources simultaneously is that an 
unrealistic model fit for a companion will spoil the fitting quality of 
the target galaxy. We develop an empirical set of criteria for iden- 
tifying bad companions such that any companion being fitted with 
(Tago/f'SO > 0.3 and n > 8.0 is identified as a bad companion, 
where aso is the semi-major axis aso = rso/ y^&T^, b/a is the axis 
ratio, and a^.^^ is the internal GALFIT error of aso. Such bad com- 
panions are usually small and/or low surface brightness galaxies 
and tend to be fit with an overestimated size, Sersic index and mag- 
nitude, resulting in a severe underestimate of the same quantities for 
the target galaxies. We examine SExtractor magnitude differences 
between bad companions (mac) and target galaxies (mT) and find 
that bad companions, as defined here, are mostly fainter than their 
corresponding targets. Actually, the majority of bad companions 
have Am = m_BC — mr > 2.5 and the distribution of Am is 
peaked at Am ~ 3.5. Thus we simply mask out all Am > 2.5 bad 
companions and refit the image. We iterate this procedure until all 
the targets have no more bad companions. About 8 percent of our 
CEN sample start with bad companions and we correct all of them 
as described above. Some of the target CENs have Am < 2.5 bad 
companions within a few arcseconds and any masking might ad- 
versely affect the fit. Hence instead we opt to exclude those target 
galaxies from our CEN sample. The fraction of CENs having such 
bad companions depends on halo mass: from 7 percent at high- 
masses to 3 percent at low-masses. In total, 32 CENs are excluded 
from our sample owing to problems with bad companions. 



3.4 Sky Uncertainty and Parameter Covariance 

Although the SDSS global sky works as a good choice of sky back- 
ground for fitting isolated galaxies, as the simulations show, the real 



sky is difficult to measure especially for CENs in dense environ- 
ments. In this subsection, we estimate the uncertainty in the SDSS 
global sky values and evaluate how this uncertainty translates into 
fit parameter errors. 

We use the difference between the SDSS global sky and mean 
sky to characterise the uncertainty in the sky measurement. SDSS 
measures the global sky as the median data counts ( ADUs) from ev- 
ery pixel of the source-subtracted frame after sigma-clipping. Be- 
sides the global sky, SDSS also provides a mean sky for each frame 
and a local sky for each detected object in the frame. SDSS mea- 
sures the background in sub-frame boxes of of 256 x 256 pixels 
centred every 128 pixels. The local sky of each object in a SDSS 
frame is an interpolation of the sub-frame background values at 
the position of the object centre. The mean of all sub-frames is the 
mean sky for a frame. 

In the left panel of Figure |4l we plot the difference between 
the SDSS global sky and mean sky for our CEN sample. The mean 
sky is somewhat larger than the global sky on average because indi- 
vidual sub-frames may suffer contamination from large and bright 
sources. We also plot the difference between the SDSS global and 
local sky measurements for our sample in the left panel. These dif- 
ferences are much larger than that between the global and mean sky 
because the local sky is heavily contaminated by large and bright 
sources. Furthermore, in the middle panel of Figure |4l we show 
that the overestimates of the local sky are halo-mass dependent for 
CENs. The local sky at the centre of massive haloes tends to be 
higher than that in less massive haloes. Since the real sky back- 
ground should be independent of the properties of haloes, this de- 
pendence implies that local sky measures for CENs suffer from in- 
creased contamination from the higher density of galaxies in larger 
haloes. In contrast, SATs display little dependence on halo mass, 
even though their local sky is also overestimated, as shown in the 
right panel of Figure|4] We attribute this lack of halo dependence to 
the fact that SATs are found over a range of halo-centric positions 
and hence at various galaxy densities. It is also possible that intra 
cluster light (ICL) at the centre of massive haloes can contribute to 
the difference between CEN and SAT sky estimates as we show in 
Figure |4] Given the large overestimates and the halo dependence 
of the SDSS local sky, we will not use it as a measurement of the 
sky background. Moreover, we conclude that the SDSS mean sky 
not only provides a robust alternative measurement of the sky back- 
ground but it also tells us in which direction the sky measurement 
might be biased. 

To study the effect of sky uncertainties on the Sersic fitting of 
actual galaxies, we first select a representative subset of 45 CENs 
from our sample that span a 3 x 3 matrix in Mhaio-w space as 
follows: 

• for log(Mhaio/h"^Mo)=[I2.0:12.5], n=[1.8:2.2], [3.3:3.7], 
[4.8:5.2] 

• for log(Mhaio/h"^M0)=[I3.O:13.5], n=[2.8:3.2], [4.3:4.7], 
[6.3:6.7] 

• for log(Mhaio/h~^M0)=[14.O:14.5], «=[3.8:4.2], [4.8:5.2], 
[6.3:6.7] 

Next, we randomly select five galaxies from each bin of the above 
matrix. For each galaxy, we Monte-Carlo sample 50 values from 
the distribution of the mean-global difference shown by the red line 
in the left panel of Figure |4] These values are added to the global 
sky and the galaxy is refit using these new background levels. This 
procedure provides a distribution of fitting parameters caused by 
the uncertainty in measuring the sky. For this subset of 45 CENs, 
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we plot distributions around the mean value of three important fit 
parameters (mtot, n and aso) in FigurejS] 

The distributions in Figure[5]clearly demonstrate a covariance 
between best-fitting parameter values and the choice of sky. The 
boxes represent our original results using the SDSS global sky. 
The series of (50) crosses, which form a short and nearly straight 
line through each box, represent the fitting results using the Monte 
Carlo sampling of the sky background as described above. The cor- 
relations between mtot and aso (bottom) and nitot and n (top) show 
the strong degeneracy of these parameters in the Sersic model. If we 
increase the sky background, the best-fitting flux decreases to keep 
the total (galaxy + sky) flux constant, while the best-fitting n and 
aso decrease. Conversely, decreasing the sky background level re- 
sults in larger best-fitting values for mtot, n, and aso. We find that 
the covariance is more severe for galaxies with higher n. This effect 
shows that the sky estimate is very crucial for producing accurate 
Sersic profile fits. This is especially true for galaxies with n > 4 
because the flat, extended wings of the profile are sensitive to sky 
uncertainties. The strength and direction of these covariances must 
be accounted for when analysing the size-luminosity relation, as we 
will discuss in ^4.2\ 



4 THE STRUCTURE OF CENS 

In this section we explore the r-band structural properties of our 
CEN sample. CENs are the most massive members of the SDSS 
groups and clusters, and are presumed to be located at the dynami- 
cal centre of the host halo. We focus on the structural shape (char- 
acterised by the Sersic index n) and the size (characterised by the 
half-light radius rso) of the Sersic profile of CENs. We study the 
relationships between these parameters and galaxy stellar masses 
and host halo masses to investigate which factor is more related to 
the structure of CENs. We also compare the structural parameters 
of CENs and SATs using our two matched SAT samples (see ij2| 
for details) and study whether the central halo location impacts the 
shape or size of a galaxy. 

To evaluate the accuracy of our measurements of galaxy struc- 
tural parameters and total flux using our GALFIT pipeline, we 
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compare our fit results with the one-dimensional ( ID) Sersic fit- 
ting p arameters based on the NYU-VAGC analysis telanton et alj 
l200i hereafter 'NYU-VAGC'). The NYU-VAGC Sersic parame- 
ters have been widely used in the study o f gala xy morphology an d 
size (e.g. lShen et al.ll2003l : lBlanton et al.ll2005l : lMaller et alfcoosh . 
We find that for high Sersic index (our GALFIT n > 3.0) galaxies, 
NYU-VAGC underestimates the Sersic index by about 1.3 and un- 
derestimates the total magnitude by about 0.4 mag compared to our 
results. There are two reasons for these underestimates: (1) NYU- 
VAGC's ID profile fitting systematically underestimates the Sersic 
in dex and total magnitu de for high n galaxies, as shown by Figure 9 
of lBlanton etal] ( |2005|) : (2) NYU-VAGC uses the SDSS local sky, 
which overestimates the background sky level and hence results in 
underestimates of n and the total magnitude. Furthermore, for low 
Sersic index (our GALFIT n < 3.0) galaxies, NYU-VAGC over- 
estimates the Sersic index but still underestimates the total magni- 
tude by about 0.2 mag. Besides using the local sky, NYU-VAGC's 
azimuthally-averaged ID fitting procedure systematically overes- 
timates the Sersic index for highly inclined galaxies, which are 
mostly disks. We also compare both our GALFIT results and the 
NYU-VAGC results with the Petrosian quantities from the SDSS 
photom etric pipeline. By following the formalism of Gra ham et al.l 
(l2005l) . we find that under our working assumptions (a Sersic model 
and a global sky), our GALFIT fitting pipeline more accurately 
measures the structural properties of SDSS galaxies then NYU- 
VAGC's. The details of this comparison are in Appendix IaI 



4.1 The Sersic Index of CENs 

With reliable and accurate measurements of galaxy structural prop- 
erties from our fitting pipeline in hand, we turn to the study of 
how the Sersic index of CENs is related to their stellar mass and 
their environment. The Sersic index is widely used to characterise 
the p rofile and concentration of galaxies in both obs ervational 
1996; 



20031: 



Graham etal. 1996; Blanton et al. 2005; MacArt hur et al.1 
de Jong et al. 2004; Allen et al. 2006) and numerical (e.g.. 



Naab & Truiilla2006. ; ,Aceves et al..2006i) studies. Some argue that 



the morphology-density relation implies that the structure of galax- 
ies is affected by the environment. Others argue that the structure 
of galaxies depends str ongly on stellar mass but only weakly on 
the environment (e.g., iHogg et al.l |2004 iKauffmann et al.l |2004 
Ivan der Welll2008l) . In this section, we explore the dependence of 
the Sersic -index distribution of CEN galaxies on both Mstar and 

Mhalo. 



4.1.1 Dependence on Stellar Mass 

In the upper left panel of Figure[6] we plot the Sersic index of CENs 
as a function of their stellar mass. The Mstar v alues for oiir CEj SI 
sample are calculated following the formula from lBell et aljjioosh : 

log[Mstar/MQ] = -0.306+1. 097"-''(g-r)-0.15-0.4(°-°Mr-4.67), (3) 

where the constant 0.15 corrects to a iKroupj i200lh IMF and 
° °{g - r) is the Petrosian col our from the NYU-V AGC shifted to 
the 2 = rest-frame using the lBlanton et alj i2003t) A'-corrections 
and co rrecting for Milky Way extinction using the lSchlegel et al.l 
( 11998 ') dust maps. The absolute r-band magnitude ^'^ Mr is ex- 
tinction and K-corrected in the same manner, but we use the to- 
tal flux from our fitting for this quantity rather than the Petrosian 
magnitude in the SDSS pipeline. We show in Appendix |A] that 
our fitting procedure better recovers the total magnitude of galax- 
ies than the Petrosian photometry. The red symbols with error- 



bars in Figure [6] show the median and the 1st and 3rd quartile 
of the n distribution in Mstar bins with a 0.25 dex width. From 
the plot we see that the median Sersic index is a strong function 
of stellar mass: low-mass CENs have low n (disk-like) and high- 
mass CENs have high n (spheroid-like). Yet, there is large scatter 
in the relation between n and Mstar , especially for CENs with 
3 X 10"' < Mstar /h^^M© < lO", where the Sersic index ranges 
from 1 < n < 6 or larger. Higher-mass CENs tend to have n > 3. 

We also show how the uncertainty and covariance of the GAL- 
FIT n values affect the n-Mstar relation in the upper left panel of 
Figure |6] (green symbols). The red symbols in the plot only take 
into account of the scatter from our fits. However, as we showed in 
i]3.4l there is a strong covariance between the fit parameters n and 
mtot and the sky uncertainty, and we use mtot to estimate stellar 
mass. To evaluate how this covariance may change the above n- 
Mstar relation, we re-calculate the relation using the results from 
the representative subset of 45 CENs in i]3.4| For each CEN in our 
sample we randomly select one CEN from the nearest five in n 
space among the representative subset, and we assume that each 
CEN has the same n and mtot covariance as its matched com- 
panion. In this way, we construct a probability distribution of n 
and Mstar (from mtot) for all 91 1 CENs. The median, 1st and 3rd 
quartiles of this distribution are shown by the green symbols and er- 
rorbars in the upper left panel of Figure (6) The new relation shows 
little difference compared with the original because the parameter 
uncertainties owing to sky are smaller than the intrinsic scatter in 
our sample (red symbols). For example, the average relative scatter 
between the 1st and 3rd quartiles of n owing to the background sky 
level uncertainty is only /S.n/n ~ 0.2, while the measured scatter 
is An/n ~ 0.75 and hence dominates. 



4.1.2 The Dependence on Host Halo Mass 

In the upper right panel of Figure[6l we show the dependence of the 
Sersic index of CENs on their host halo mass. Recall that Mhaio is 
calculated by matching the rank of the total stellar mass of groups 
and clusters with that of dark matter haloes from numerical simu- 
lations (see Y07 for details). We find that n depends only weakly 
on Mhaio and that the scatter of the n-Mhaio relation is large. For 
CENs in haloes with log[Mhaio/h"^ (Mq)] < 12.5, we find a me- 
dian value of n = 3 and a relative scatter of An/n ~ 1.3. For 
massive haloes with log[Mhaio/h~'^(M0)] > 14.0 the median is 
71 = 5 and the relative scatter is An/n ~ 0.6. 

The 71-Mhaio relation, although weak, suggests that the halo 
mass may also affect the structure of CENs. However, from the 
right panel of Figure[T]we know that the stellar mass of CENs also 
depends on the halo mass in the sense that CENs in massive haloes 
tend to have larger stellar masses than CENs in smaller haloes. 
Given this dependence and the strength of the n-Mstar relation, 
it is tempting to rule out any dependence of n on halo mass. To 
address this, we attempt to remove any Mstar -Mhaio dependence 
from both the n-Mstar and n-Mhaio relations. 

First, in the lower left panel of Figure |6l we plot the n-Mstar 
relations for CENs in five halo mass bins, each 0.5 dex in width. All 
the bins with less than six CENs are excluded to get better statis- 
tics. The roughly Mhaio -independent relations all have a slope and 
amplitude that is similar to the single n-Mstar relation for the full 
sample (red symbols in the upper left panel). There is some evi- 
dence that the relations for different Mhaio bins are somewhat off- 
set along the Mstar direction in the sense that CENs in less massive 
haloes tend to have larger n than their counterparts in more massive 
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Figure 6. Upper left: the relation between Sersic index and stellar mass for our CEN sample. Black crosses represent the best-fitting n for individual CEN 
galaxies. The red line connects the median (open squares) of the n distribution in each Mstar bin. The green Hne and open circles show the median once the 
n and mtot fitting uncertainties owing to sky are folded into the distiibution (see text for details). The eiTorbars show the first and third quartiles. Upper right: 
the relation between the best-fitting Sersic index and halo mass for our CEN sample. Black crosses and red lines with errorbars are done as in the left panel. 
Lower left: the n-Mstar relation split into different bins of fixed halo mass as indicated by key. Lower right: the n-Mhaio relation split into different bins of 
fixed stellar mass as indicated by the key. In both lower panels, the symbols and errorbars provide the median, 1st, and 3rd quartiles. In the lower panels, only 
the best-fitting results are used. 



haloes. But since the scatter in the relations are large, we cannot 
draw any firm conclusions. 

Similarly, we study the n-Mhaio relations of CENs in five dif- 
ferent Mstar bins, each with a width of 0.4 dex, as shown in the 
lower right panel of Figure |6] We find that the n-Mhaio relations 
for each bin of fixed galaxy mass is different from the single n- 



Mhaio relation for the full sample. For the most part, the relations 
are all flat, i.e. each CEN of a given Mstar has a constant n, within 
the scatter, independent of its host halo mass. We note that several 
of the fixed stellar mass relations have a small negative slope when 
considering only the median values. This slight trend is a manifes- 
tation of the small offsets among the n-Mstar relations in different 
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Mhaio bins in the lower left panel. A much larger sample is re- 
quired to validate whether the profile shape of CENs, as measured 
with Sersic fitting, has a small second-order dependence on halo 
mass. 

4.1.3 Discussion 

We study the distribution of n for our CEN sample and its depen- 
dence on Mstar and Mhaio and find that the Sersic profile shape 
of CENs is strongly correlated with Mstar but only weakly (or 
not at all) correlated with Mhaio- Low-mass CENs tend to have 
shallower, disk-like profiles, while massive CENs have steeper, 
more spheroid-like shapes. This n-Mstar relation holds for CENs 
from different Mhaio bins with almost the same slope and ampli- 
tude. On the other hand, CENs have n values that depend very 
weakly, if at all, on global environment as defined by Mhaio, 
which is consisten t with other observations ( Kauffmann et al . 20041 : 
Ivan der We]||2008h . This correlation disappears (or even becomes 
an anti-correlation) if we divide our CEN sample into different 
Mstar bins. This suggests that Mhaio has almost no effect on the 
shape or concentration of CENs. The key factor that determines n 
of a CEN is its stellar mass. Given the relationship between the 
masses of CENs and their host halo masses, the phenomenon of 
CENs in massive haloes having a more concentrated structure (as 
implied by morphology-density relation) can be explained simply 
by (1) an intrinsic n-Mstar correlation, and (2) massive CENs liv- 
ing preferentially in massive haloes (i.e., mass segregation). 

We find low-n (disk-like) CENs in haloes spanning a large 
range in halo mass (from 10^^ to lO^*h~^M0), which rules out 
a distinct (special) Mhaio for producing spheroids. Likewise, find- 
ing high-n (early-type) CENs over the same range of environments 
suggests that major mergers, or some other processes that trans- 
forms disk galaxies into spheroid-like galaxies, occur at the centres 
of haloes with a wide range of masses. Finally, it is unclear whether 
the halo plays a small, secondary role in determining the overall 
profile shape of its CEN as hinted at in the bottom panels of Fig- 
ure[6] The Mstar -Mhaio relation in Figure 1 shows that CENs of a 
fixed stellar mass reside in haloes with a range of masses. We spec- 
ulate that the tendency for CENs of a fixed stellar mass to be more 
concentrated in smaller haloes and less so in larger haloes may tell 
us something about the recent accretion history of the host and its 
impact on the stellar mass growth of the CEN. Halos grow through 
the accretion of, and the merging with, other haloes. Briefly, con- 
sider two haloes (A and B) that were equal in mass and contained 
similar central disk galaxies in the recent past. Since that time both 
haloes have doubled in mass, halo A by a single major merger with 
a comparable halo, and B by accreting many minor subhaloes. Ow- 
ing to dynamical friction timescales, the disk CEN in halo A will 
both (i) grow in mass faster than halo B's CEN and (ii) be trans- 
formed into a spheroid because it is doomed to experience a major 
merger with the CEN of the merging halo. Conversely, the mass of 
the CEN in halo B will take much longer to grow by the occasional 
minor merger with small infalling companions, and this process is 
much less likely to destroy its disk morphology. Therefore, when 
comparing haloes of the same mass, galaxies that grew through 
major mergers will be more massive and more concentrated. Like- 
wise, when comparing CENs of the same stellar mass, those which 
grew through major mergers (i.e., more concentrated) will reside 
in lower mass haloes on average more than their disk-dominated 
counterparts as we see in the lower left panel of Figure [6] Numer- 
ical modelling and a more-detailed analysis on a much larger data 
set could test these predictions. 
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Figure 7. The distribution of Sersic index for our visual inspected early- 
(solid line) and late-type (dotted line) CENs. The Sersic index comes from 
our best fits. 



4.2 The Size of CENs 

In addition to the Sersic index, the half-light size is an important 
characteristic of galaxy structure. It is well-established that galax- 
ies follow well-defined size-luminosity (rso-i) and size-stellar 
mass (rsQ-Mr) rela tions (e.g. Shen et al. 2003; Bernardi et al. 200"^; 
iDutton et alj|2007h . which are commonly used to constrain galaxy 
formation and evolution theories. In this subsection, we study the 
size-luminosity relation and the size-stellar mass relations of our 
CEN sample and compare our results with those obtained by oth- 
ers. Our CEN sample contains both early and late-type galaxies. 
Thus, we visually inspect each galaxy and divide the sample into 
two types based on whether spiral disk structure is present (late- 
type) or not (early-type). In Figure|7|we plot the distribution of our 
best-fitting Sersic index for our visual inspected early- (solid line) 
and late-type (dotted line) CENs. Unlike employing a sharp cut at 
n — 2.5 as in many studies, the majority of our early-type CENs 
have n > 3.5, while late-type CENs have n < 3.5. In what fol- 
lows, we will discuss the early and late-type size-luminosity and 
size-stellar mass relations separately. 

In FigurejS] we show the size-luminosity relation of early- (the 
upper left panel) and late-type (the lower left panel) CEN galax- 
ies. Here we use the half-light radius rso (in units of h~^ kpc) 
to represent the size of a galaxy. We have rso = x x 
yjb/a X dA{z), where aso and b/a is the semi-major axis and 
the axis ratio of our best-fitting model of the galaxy, dA{z) is 
the angular diameter distance at redshift z of the galaxy. We cal- 
culate the absolute magnitude using the total m agnitude from our 
Sersic fit, K-l-E corrected to 2; = 0.1. Following I Shen et alj ( l2003l 
hereafter S03), who also measured the tsq-L or rso-Mr relation 
for SDSS galaxies, we use linear regression to fit the r^o-L rela- 
tion, i.e. Iog[r5o/(h~^kpc)] — — 0.4aMr + /9. We find slopes of 
Q = 1.02 ± 0.03 for early-type CENs and a = 0.88 ± 0.04 for 
late-type CENs. The linear fits are shown by thick solid lines in 
each panel. 

At the bright end, Mr — 5 log h < —22, the relation steepens. 
The origin of this steepening is unclear. It could be caused by the 
covariance between the semi-major axis aso and the total magni- 
tude mtot, as shown in the upper panel of Figure [5] Our tests in 
i|3.4| show that there is a very strong covariance between aso and 
mtot , owing to uncertainties in the measurement of background sky 
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Figure 8. The size-luminosity (left) and size-stellar mass (right) relations for early-type (upper) and late-type (lower) galaxies in our CEN sample. The 
thick solid line in each panel is the best fit linear relation for each morphological type over the whole luminosity/stellar mass range. The dashed line in each 
panel shows the linear relation by fitting the faint end (Mr — 51og(h) > —22). For comparison, the thin dotted line in each panel is the relation from S03 
(see the text for details). The arrow in the upper right panel shows the direction and amplitude of the average covariance between Mr and rso for bright 
(Mr — 51og(h) < —22) CENs owing to background sky uncertainties (see i|3.4| for details). 



level, and that this covariance becomes stronger for bright and large 
galaxies. The covariance between aso and mtot will produce a co- 
variance between 7-50 and Mr in Figure [8] since Mr is calculated 
from mtot- The direction of the average covariance at the bright 
magnitude end, as shown by the arrow in the upper left panel of Fig- 
ure[8l is almost parallel to the slope of the relation at the bright end. 
Moreover, the amplitude of the covariance increases for brighter 
and larger galaxies as seen in Figure [5] Hence the covariance be- 
tween size and total magnitude in the profile fitting could contribute 
to the slope steepening at the bright end. 

To account for any bias at the bright end, we refit exclud- 
ing sources with Mr — 51ogh < —22 and find slopes of a = 
0.82 ± 0.06 (early-type) and a = 0.57 ± 0.04 (late-type). For 
reference, a fit to the r^o-L relation of early-type CENs brighter 
than -22 has a slope of a = 1.32 ± 0.08. We don't fit the bright end 
of late-type CENs, because there are only several galaxies. Using 
the fit over the faint end (Mr — 5 log h > —22), we find an early- 



type rsQ-L relation with a much steeper slope than the a — 0.65 
of 803 (the thin dotted line in the upper left panel of Figure [Sj. 
Our fit to the r^o-L relation for late- type CENs fainter than -22 is 
also qualitatively steeper than that of S03 (comparing the dashed 
and dotted line in the lower left panel), who fit the r^o-L rela- 
tion with a four-parameter model. We note that the definition of 
early/late types in S03 is different from ours. They classified galax- 
ies with a Sersic index n > 2.5 as early- type and the others as 
late-type. The different definitions of early/late-type between S03 
and us could contribute to the measured slope differences. We also 
try the same early/late-type classification as S03 by using Sersic 
index from our best fits and refit the r^o-L relation for the faint 
end CENs. We find a — 0.78 ± 0.05 for the slope of early-type 
CENs, which becomes flatter due to containing more low Sersic 
index galaxies than our visually inspected early-type sample but is 
still steeper than the slope of S03. There are two other possible rea- 
sons for the large discrepancies in slope between the r^o-L relation 
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of S03 and our result: (1) S03 use the results of an early version 
of the NYU-VAGC fitting, which underestimates both the total flux 
and the size of galaxies, as shown in Appendix IaI (2) S03 covers 
a wide range of luminosity and hence contains many faint (fainter 
than the faintest one in our sample) galaxies. The slope of the rso- 
L relation could have a smooth transition from a flat one at the 
faint gala xy end to a steep one at the luminous galaxy end, as sug- 
gested bv lPesroches et all 12007). Including many faint galaxies in 
the S03's sample can flatten the average slope of the whole sample, 
making their slope smaller than ours. 

We also compare our results with other studies of the rso- 
L relation of BCGs, a subse t of massive early-ty pe galaxies at 
the centres of large clusters. Bemardi et alj ( 120071) fit a de Vau- 
couleur s model to SPSS i mages of BCGs from the C4 cluster cat- 
alogue jMiller et alj|2005l) and study the r^o-L relation by using 
the half-light radius of their best-fitting models. They find a slope 
of a = 0.89 for early-type BCGs over the luminosity range of 
Mr - 51og(h) = [-20.3 : -23.3]. This slope is a bit shallower 
than our slope (~ 1.0) for early-type CENs over the same lumi- 
nosity range. The difference could be attributed to the fact that 
they fit BCGs with a de Vaucouleurs profile rather than a Sersic 
profile, which yields a smaller size a nd hence results in a flatter 
rso-L relation, as shown in their paper. lvon der Linden et al.l j2007l) 
performed isopho tal photometry on the SDSS BCGs from the C4 
cluster catalogue jMiller et al1l2005l) and found a = 0.65 ± 0.02 
for BCGs over the luminosity range Mr — 51og(h) — [—20.3 : 
—23.3]. However, the isophotal photometry technique could miss 
the extended outer parts of galaxy l ight profiles and h ence result 
in underestimates of both rso and L. lLauer et al. II2OO7I) combined 
surface photometry presented in several HST imaging programs for 
219 early type galaxies and fit them using a de Vaucouleurs pro- 
file. They found a = 1.18 ± 0.06 for galaxies with Mv < -21 
and that the r^o-L relation changed from a flat slope to a steep 
slope at around My = —22. This slope is close to our slope for 
CENs brighter than Mr — 51og(h) — —22. Also, their slope be- 
haves qualitatively like ours, with a stee p slope at the bright e nd 
and a more shallow slope at the faint end. lGonzalez et al] j2005i) fit 
Sersic profiles to the I-band images of 24 luminous BCGs residing 
in large clusters and found a = 1.8 ± 0.2, but these fits included 
significant ICL. During our tests in WiA\ we showed the sensitivity 
to the background sky estimate and the covariance of the size and 
magnitude from profile fitting. This could become more acute for 
high luminosity CENs, where distinguishing the ICL from the ex- 
tended outer light bound to a galaxy becomes increasingly difficult. 

Overall, although it is well known that size strongly correlates 
with luminosity, the slope of this relation is hard to determine pre- 
cisely. As we summarised in Table [T] authors with different sam- 
ples, measurement techniques, definitions of size and/or even the 
criteria to split their samples into early and late types obtain dif- 
ferent slopes. What's more, a small uncertainty in background sky 
translates into a strong covariance between size and luminosity and 
hence can change the slope, especially for bright and large galaxies. 
Based on our analysis of the robustness and uncertainty of our fit- 
ting, and combined with previous work, we propose slopes for the 
rso-I/ relation of q ~ 0.9 and a ~ 0.6 for early-type and late- type 
CENs, respectively. 

Closely related to the size-luminosity relation of CEN galax- 
ies is their size-stellar mass rso-Mstar relation. For comparison, 
we plot the rso-Mstar relations for our early- type (upper right) and 
late-type (lower left) CENs in Figure [S] We find slopes of a = 
0.90 ± 0.02 (early-type) and a = 0.47 ± 0.03 (late-type) by using 
the relation log[r-5o/(h"^kpc)] = alog[Mstar/(h"^M0)] + to 



fit the whole luminosity/stellar mass range of our CEN sam- 
ple. Similar to the r^o-L relations, our rso-Mstar relations are 
steepened by the covariance between size and total magnitude 
in the profile fitting. We also refit the relations with bright 
(Mr - 51og(h) < -22) CENs excluded and find slopes of a = 
0.68 ± 0.03 (early-type) and a = 0.35 ± 0.03 (late-type). Our 
slopes of faint end (Mr — 51og(h) > —22) CENs are steeper than 
those in S03 (thin dotted lines in each panel). As discussed above, 
differences in the fitting procedures could be responsible for the 
difference. Also we note that the rso-Mstar relation in S03 is de- 
rived by using the z-band rso ■ 

4.3 Comparison Between CENs and SATs 

One of our two primary questions regarding the environmental de- 
pendence of morphological transformation is whether or not the 
group centre is a special place for determining the structure of 
galaxies. In this subsection, we compare the structural properties 
of CEN galaxies with those of galaxies selected from the two SAT 
samples, which are comparable to our CEN sample. First, we con- 
sider SATs matched in stellar mass with our CEN sample, the SAT 
sample SI. Second, we study SATs matched both in stellar mass 
and colour, the SAT sample S2. (see Sj2]for details). We address our 
above question by testing if CENs and SATs are two distinct popu- 
lations of galaxies possessing different structural properties. SATs 
were CENs before being accreted by a larger halo, thus, differences 
between these galaxies probe the impact of local environment on 
SAT-specific transformation processes. 

Many studies have demonstrated structural differences be- 
tween brightest cluster galaxies (BCGs) and non- BCGs. For ex- 
ampl e, BCGs are found to have larger sizes (e. g., iBemardi et al.l 
i2007i : lvon der Linden et alJl2007l,lLiu et al.ll2008h and steeper light 



profiles (e.g., lGraham et aljl996h compared to other massive early 
type galaxies (ETGs) in the same cluster. These results are inter- 
preted as differences between CEN and SAT galaxies and, as such, 
are thought to indicate unique formation histories. However, it is 
important to keep in mind that the BCGs in these studies represent 
a special subset of all CENs; namely they are the most-luminous 
and highest-mass galaxies found at the centres of massive clusters 
and are the tip of the galaxy group mass function. Moreover, a com- 
parison of BCGs with other morphologically similar cluster mem- 
bers is effectively a comparison of galaxies of different masses. In 
general, the structural properties of ETGs have a s mooth transition 
as their luminosity and/or stellar mass changes jPesroches et al.l 
12007). Therefore, it is unclear whether differences between BCGs 
and non-BCGs are intrinsic and originate from separate formation 
mechanisms or are simply a reflection of their different masses. To 
avoid this selection effect, we compare the structure of CEN and 
SAT galaxies that are matched in stellar mass. 

4.3.1 Sersic Index n 

We compare the n-Mstar relations of CENs and SATs in Figure 
|9] The relation for CENs is reproduced from Figure |6l and the n- 
Mstar relations of SAT sample S 1 and S2 are measured in an iden- 
tical fashion as for the CENS, i.e. their stellar masses are based 
on their best-fitting total magnitudes from our Sersic fitting analy- 
sis and the error bar represents the first and third quartile in each 
Mstar bin. We also repeat our self-matching test in Sec. 12.21 and 
find that the scatter of Sersic index medians due to our match- 
ing scheme is much smaller than the sample scatter So we de- 
cide to choose the scatter of the sample (as shown by the first 
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Table 1. Slope of the r^o-L relation for early-type galaxies from various authors 



Work Sample a Aa Luminosity range Model 



This work 


Y07 CENs 


1.02 


±0.03 


Mr 


- 51og(h) = [-19 


: -24] 


Sersic 


This work 


Y07 CENs 


0.82 


±0.06 




Mr - 51og(h) > - 


22 


Sersic 


This work 


Y07 CENs 


1.32 


±0.08 




Mr - 51og(h) < - 


22 


Sersic 


Shen et al. (2003) 


SDSS galaxies 


0.65 




Mr - 


51og(h) = [-15.3 


: -23.3] 


Sersic 


von der Linden et al. (2008) 


C4 BCGs 


0.65 


±0.02 


Mr - 


51og{h) = [-20.3 


: -23.3] 


isophotal photometry 


Bernardi et al. (2007) 


C4 BCGs 


0.89 




Mr - 


51og(h) = [-20.3 


: -23.8] 


de Vaucouleurs 


Lauer et al. (2006) 


219 galaxies 


1.18 


±0.06 




Mv < -21 




de Vaucouleurs 


Gonzalez et al. (2006) 


24 BCGs 


1.8 


±0.2 




Mi < -24 




Sersic 
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Figure 9. Comparison of the n-Mstar relation for CEN galaxies with that 
of SATs matched in stellar mass (SAT sample S 1 , upper panel), and with 
SAT galaxies matched to our CEN sample in both colour and stellar mass 
(SAT sample 32, lower panel). In each panel, solid lines with filled circles 
represent the median values of the CEN relation and dot-dashed lines with 
open squares show the SAT samples. The en'orbars give the first and third 
quartiles of the n distribution in each Mstar bin as in Figure |6] and bins 
with less than six galaxies are excluded. 



and third quartiles in the plot) to represent the confidence inter- 
val whenever comparing medians of measurements from two sam- 
ples. Note that we match our CEN and SAT samples using the stel- 
lar mass estimated from the Petrosian magnitude in ^ but plot 
the result using the stellar mass estimated from the total magni- 
tude in our Sersic fits. The only possible concern about this pro- 
cedure is that it may change the colour difference between CENs 
and SATs in the SAT sample S2 at given stellar mass (see the 
lower right panel of Fig. To evaluate the magnitude of this 
effect we calculate the difference between the median colour of 
CENs and SATs in the SAT sample S2 in the same stellar mass 
bin with a width of 0.5 dex from log[Metar/(h-2M0)] = 9.5 tO 
log[Metar/(h-2M0)] = 12.0, where now the stellar mass is the 
one computed by using our Sersic total magnitude. We find a very 
small difference: SATs are redder by 0.015. However, this differ- 
ence is still well within the measurement uncertainties of colour 
(±0.03) and hence our claim that the SAT sample S2 is matched 
with the CEN sample in both Mstar and colour is still valid. As for 
the SAT sample SI, as long as we compare CENs and SATs in same 
stellar mass bin, the modification of stellar mass does not change 
the results. 

In the lower panel of Figure |9] we see that the n-Mstar re- 
lation of SATs is almost identical to that of CENs matched in 
BOTH stellar mass AND colour. The only exception is for SATs 
in our highest-mass bin, which have a somewhat higher median 
Sersic value than CENs of similar mass and colour. We note that 
this bin has the smallest number of SATs (13), so small number 
statistics may account for the difference in the medians. How- 
ever, if we release the constraint on colour matching, low mass 
SATs (log[Mstar/(h"^MQ)] < 10.75) tend to have a somewhat 
higher median Sersic index than their CEN counterparts (upper 
panel), although the discrepancy is within the scatter. In contrast, 
the more massive SATs have similar n-Mstar relations as CENs 
matched only in stellar mass, except again for the highest-mass 
SATs. Our results are in g ood agreement with those obtained by 
Ivan den Bosch et al] ( |2008|) . who used concentration defined by 
the ratio of SDSS radii containing 90 percent and 50 percent of 
the Petrosian flux, rather than n, to describe the profile shape of 
galaxies and found that low mass SATs are redder and slightly 
more concentrated if match ed with similar mass CENs. In addition, 
Ivan den Bosch et ai]( l2008t) found that the concentration difference 
goes to zero when matched in colour and M as we find here. 



4.3.2 Size 

Besides the n-Mstar relations, we use the rso-L and rso-Mstar re- 
lations to compare the structural properties of CEN and SAT galax- 
ies. We plot the sizes, luminosities, and stellar masses of our CEN 
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Figure 10. Comparison of CEN (red points) and SAT (blue points) galaxy size-luminosity (left column) and size-stellar mass (right column) scaling relations. 
These SATs (sample SI) have similar stellar masses as our CEN sample. The black (green) lines and errorbars show the median, first and third quartiles of the 
CEN (SAT) distribution in each luminosity or stellar mass bin. Both CENs and SATs are separated into early-type (upper panels) and late-type (lower panels) 
galaxies (see text for detail). 



sample from Figure [8] as red points in Figures [TO] and Figure [TT] 
The blue points represent the data for our SAT samples Si (Fig. 
llOl l and S2 (Fig.lilt. The size, absolute magnitude and stellar mass 
of SATs are calculated in the same way as CENs as described in 
^431 and ^42l Also, as in ^431 we divide the SAT samples SI and 
S2 into early-type and late-type galaxies based on our visual in- 
spection. However, instead of computing the slopes of the CEN and 
SAT scaling relations, we compare the half-light size distributions 
(median, first, and third quartiles) of each sample in narrow bins of 
luminosity (0.25 mag wide for the r^o-L relation) and stellar mass 
(0.2 dex wide for rso-Mstar )• In this manner, we directly compare 
CENs with SATs from each sample only in regions of overlapping 
luminosity or stellar mass. Hence, we avoid comparing scaling re- 
lation slopes based upon samples that span different luminosity and 
stellar mass ranges. This helps to avoid the effects of small number 
statistics at the bright/massive and the faint/low-mass ends, which 
can bias the scaling relation slopes. 

The comparisons in Figures[TO]and[TT]show that for early type 
galaxies, whether or not they are matched in colour, CENs and 
SATs display almost no difference in their size distributions. The 
only differences occur at the bright/massive ends of the SAT rela- 



tions, where both SAT sample Si and S2 suffer from small num- 
ber statistics. For late-type SAT galaxies, the two samples also have 
size-luminosity relations that are similar to the CENs (the lower left 
panel in both figures). SATs have smaller median sizes than CENs 
when we only match in stellar mass (the lower right panel of Fig- 
ure [Tol l, but this difference disappears when we compare to SATs 
matched in both colour and stellar mass (the lower right panel of 
Figure [TT). 



4.3.3 Discussion 

Using the n-Mstar, t^q-L and rso-Mstar relations, we compare 
galaxies in our CEN sample to SATs of the same stellar mass. 
We find two basic differences between CEN and SAT galaxies: (1) 
low mass (log[Mstar/(h~^M0)] < 10.75) SATs have a shghtly 
higher median Sersic index compared to CENs of the same stellar 
mass; and (2) low-mass, late-type SATs have smaller median sizes 
compared to their same-mass CEN counterparts. Our findings are 
in goo d qualitative agreement with the results of IWeinmann et all 
( |2008|) for a much larger SDSS sample. Using the NYU-VAGC 
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Figure 11. Same as Figure fTol except comparing CENs with SATs of similar stellar mass and colour. 



half-light radii and concentrations for over 10^ galaxies, they found 
that late-type SATs are smaller and more concentrated than late- 
type CENs of the same stellar mass. We note, however, that their 
size measurements suffer from a systematic underestimate owing 
to an overestimate of the background sky levels (as we discuss in 
Appendix [Aj. 

We also find that the above two basic differences disappear 
when we compare CENs and SATs of the same stellar mass and 
optical colour. Moreover, we find no structural differences between 
high-mass (log[Mstar/(h-^M0)] > 10.75) CEN and SAT galax- 
ies, which tend to have early-type (spheroid-dominated) morpholo- 
gies, or between early-type CENs and SATs, in general. 

For lower-mass (log[Mstar/(h"^M0)] < 10.75) SATs, the 
minor differences with CENs of similar stellar mass, and the lack 
of any difference with CENs of the same mass and colour, can be 
understood in terms of our selection criteria for the SAT sample S 1 
and S2. According to the hierarchical scenario of structure forma- 
tion, larger haloes housing groups and clusters are formed through 
the accretion and merging of smaller haloes. Thus, all SATs were 
once the CEN galaxy in a smaller halo. Under this assumption, low- 
mass SATs tend to be redder on average than CENs of similar stel- 
lar mass because SATs have had their star formation quenched by 
environmental processes once they became non-CEN members of a 



larger halo jvan den Bosch et aill2008l) . In addition, as Weinmann 
et al. (2008) point out, quenching will also cause a moderate in- 
crease in the concentration and decrease in the size of late-type 
SATs, as we find here when comparing similar mass SATs and 
CENs. It is important to keep in mind that the average structural 
differences are small and, as such, SAT quenching cannot be used 
to explain the major morphological transformation of many disks 
into spheroids that is required to produce the strong morpholog- 
ical bimodality between blue and red galaxies. When we restrict 
our SAT selection to match lower-mass CENs in both stellar mass 
and colour (SAT sample S2), we preferentially choose blue SATs. 
Lower-mass CENs are typically blue (Figure[TJ and their colour is 
associated with late-type morphology (disk-like) and on-going star 
formation. Therefore, we argue that blue SATs are likely exam- 
ples of newly accreted SATs, i.e. recent CENs of the most recently 
accreted subhaloes. If true, new SATs have not been members of 
larger haloes long enough to alter their colour and structural prop- 
erties. This line of reasoning agrees with the absence of structural 
differences we find for the CENs and SATs in SAT sample S2. 

When we consider higher-mass 

(log[Mstar/(h-^M0)] > 10.75) SATs (either sample SI or 
S2), they have the same red colour (see lower left panel of Figure. 
O and highly-concentrated, large Sersic index profiles associated 
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with early-type morphologies as CENs of similar mass. While 
massive SATs are rare compared to their CEN counterparts, they 
are structurally indistinguishable. In a broader sense, we find that 
all morphologically early-type SATs and CENs have the same size 
scaling relatio ns and that any reported differences between CENs 
and SATs (e.g. lOraham et al.ll996l ; lLiu et al.l20o3) are actually the 
result of comparing two populations, e.g. BCGs and non-BCGs, 
with different stellar mass distributions. Owing to the similar red 
colours of massive CENs and SATs, there is no way to discern 
recent arrivals from long-term members in larger haloes, but it is 
clear that the transformation into spheroids does not depend on 
becoming a SAT. If that were true, we would expect massive CENs 
to be disk-like when they are clearly otherwise. Rather, we argue 
that the strong morphological transformation from disk to spheroid 
occurred at an earlier time when a massive SAT was the CEN of 
a smaller halo and that the local environment had no additional 
impact on the structure of high-mass spheroids. We find that 
the difference between disk-dominated and spheroid-dominated 
structure is more directly related to the stellar mass of a galaxy. 
Clearly, there is some relationship between the mass of CENs and 
their host halo mass (Figure [T), but further study is required to 
understand whether or not any aspect of the environment plays an 
important role in the transformation of disks and the production of 
high-mass spheroids. 



5 SUMMARY 

In this paper, we study how the structural properties of central 
galaxies (CENs) in groups and clusters depend on galaxy stellar 
mass, global environment (group halo mass), and local environment 
(central/satellite position wit hin the host halo) . We select from the 
SDSS DR4 group catalogue jYang et al]|2007l) a statistically repre- 
sentative sample of 911 CENs whose host halo masses span from 
lO" to lO^^h~^M0. We use 2D Sersic model fits to quantify the 
shape (Sersic index) and size (half-light radius) of each galaxy. To 
this end, we establish a well-tested, GALFIT-based pipeline to fit 
Sersic models to SDSS imaging data in the r-band. We summarise 
our main findings below. 

We thoroughly test the performance of our GALFIT pipeline 
on simulated and real SDSS galaxy image data. Our 2D fitting re- 
covers the structural properties of simulated galaxies with no bias, 
unlike the one-dimensional fits to azimuthally-averaged data em- 
ployed for the NYU-VAGC that systematically underestimate the 
total flux, size and Sersic index of higher-n profiles. For galaxy 
profile fitting, we also demonstrate that the SDSS global sky is pre- 
ferred over the SDSS local value as a background level measure- 
ment. We compare our fitting results with those from the NYU- 
VAGC and find that our fits include light from the outer parts of 
galaxies, which is missed when an overestimate of the (local) sky 
background is used. We test how this background uncertainty trans- 
lates into a systematic uncertainty in the fitting parameters owing 
to a strong covariance between Sersic index, total magnitude, and 
half-light size. This covariance affects bright and large galaxies 
more and could contribute to the apparent steeping in the slope 
of the size-luminosity and size-stellar mass relations at the bright 
(massive) end. 

We find that the Sersic index of CENs depends strongly 
on Mstar, but weakly or not at all on Mhaio- The depen- 
dence on stellar mass is in the sense that low mass galax- 
ies (Mstar < lO^°'^h~^M0) have lower, disk-like indices (n ~ 
2.0), while massive galaxies (Mstar > lO" °h~^M0) have higher. 



spheroid-like (~ 5) indices. Over a large range in Mhaio , from 
small groups to large clusters, any change in the n distribution 
of CENs is likely the result of the correlation between Mstar and 
Mhaio- The fact that spheroidal CENs are found at all group masses, 
and the lack of a strong n dependence on Mhaio , both rule out a 
distinct halo mass for producing spheroids. Moreover, the strong 
dependence of n on Mstar suggests that Mstar is the key factor in 
determining the shape of CENs. 

Similar to the light profile shape, the half-light size of CENs 
depends on galaxy stellar mass and luminosity. We separate our 
CEN sample into early and late-type galaxies by visual inspection 
and we find a r^o-L slope of a ~ 0.83(0.62) for early-type (late- 
type) galaxies with —22 < Mr — 51og(h) < —20. We also com- 
pare our rso-i slope for early-type CENs with those from other 
studies and find that there is fairly large discrepancy. This discrep- 
ancy could result from several factors including different samples, 
size measurement techniques, or early-type galaxy definitions. 

To study whether the structural properties of CENs depend 
on their special position at the centre of the gravitational po- 
tential well, we compare their shapes and sizes with those of 
non-CEN satellite (SAT) galaxies. We find that low mass (< 
lO.O^^'^^h^^MQ) SATs have somewhat larger median Sersic in- 
dices compared with CENs of similar stellar mass. In addition, low 
mass late-type SATs are moderately smaller in size than late-type 
CENs when matched in stellar mass, but no size differences are 
found between early-type CENs and SATs. We find no structural 
differences between SATs and CENs when they are matched in both 
optical colour and stellar mass. The small differences in the sizes 
of low-mass, late-type CENs and SA Ts are consistent with SAT 
que nching as found by oth ers (e.g., in Ivan den Bosch et al.l ( 1200 Sh 
andl Weinmann et al.U 20Q8')). The similarity in the structure of mas- 
sive SATs and massive CENs demonstrates that the local environ- 
ment has no significant impact on the structure of a massive galaxy 
that enters a denser environment and that these two populations are 
morphologically indistinguishable. 

We conclude that Mstar is the most fundamental property in 
determining the basic structural shape and size of a galaxy. In con- 
trast, the lack of a significant n-Mhaio relation rules out a clear 
distinct group mass for producing spheroids. This fact, combined 
with the existence of spheroid CENs in low-mass and high-mass 
groups, suggests that the strong morphological transformation pro- 
cesses that produce spheroids must occur at the centres of groups 
spanning a wide range of masses. 
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Figure Al. Differences between the NYU-VAGC and GALFIT Sersic parameters from r-band SDSS data as a function of the GALFIT Sersic parame- 
ters for our total sample of 1657 CEN+SAT galaxies. The lines show the quartiles of each distribution as in Fig. [3] Black points are for galaxies with 
|Asky| > 0.1 ADU and red points for I Asky[ < 0.1 ADU, where Asky = local sky - global sky. Note, in the /e/? panels it is clear that the NYU-VAGC 
fitting has an n = 6 limit. 



APPENDIX A: COMPARISON WITH NYU-VAGC Sersic 
FITS 

iBlanton et all feOOSl hereafter 'NYU-VAGC') fit Sersic models to 
the azimuthally- averaged ID profiles o utput by the SDSS photo- 
metric pipeline jStoughton et alj|2002h for all SDSS DR4 galax- 
ies meeting the Main sample criteria. Tests show that the NYU- 
VAGC Sersic fitting does well for simulated galaxies with an input 
Sersic index nin < 2, a small size or a faint magnitude. But for 
riin > 2 simulations, the NYU-VAGC fitting systematically under- 
estimates n, 7-50, and total flux [see Figure 9 in lBlanton et al. ||2003) 
for details]. For example, for a simulated galaxy with ni„ = 4 
the NYU-VAGC fitting underestimates these parameters by 15 per- 
cent. For comparison, our GALFIT fitting results for 850 sim- 
ulated Sersic galaxies placed in SDSS images show very little 
bias for ni„ > 4 galaxies (see Figure [3]l- Note that we adopt 
the global sky value from the SDSS image header in our GAL- 
FIT profile fits, while NYU-VAGC uses the local sky level. It has 
been reported that the SDSS pipeline overestimates the local sky in 



dense enviro nments fcauer et al j2007l : lvon der Linden et al.ll2007l : 
lAdelmari-M cCarthv et al. 2008). As such, Sersic fits based on over- 
estimates of the sky may result in fainter magnitudes, smaller sizes 
and/or lower Sersic indices as we demonstrate in il3.4l 

Figure lAll shows comparisons between our fit results and 
those from NYU-VAGC for our total sample of 911 CENs plus 
746 SATs from our SAT sample S2 (see ij2]for details). Below, we 
discuss in detail the discrepancies between the two fits for galaxies 
with n > 3 from GALFIT, and then for those with lower Sersic 
indices. 

Al High Sersic Galaxies 

For n > 3 galaxies, the NYU-GALFIT parameter discrepancy 
for real galaxies in Figure |aT] follows a similar trend as those be- 
tween t he input and fit para meters for simulated galaxies in Fig- 
ure 9 of lBlanton etahl ( l2005h . but with an increased amplitude. For 
example, the Sersic difference (lower left panel) grows by about 
An ~ 1.3 over the interval 3 < n < 6, compared with An ~ 0.6 
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over the same interval in the iBlanton et all ( l2005h simulations. 
Here two factors are at play: one is the systematic underestimate 
of NYU-VAGC's ID fitting procedure for steep Sersic profiles as 
demonstrated in their test fitting using simulations. The second fac- 
tor is the difference between the sky levels used in each procedure. 
In Figure IA2i we attempt to separate these two factors by split- 
ting the whole sample into GALFIT Sersic index and relative sky 
difference bins, where the relative sky difference is the difference 
between the local and the global sky normalised by the Petrosian 
surface brightness of the galaxies (Ipetro = fpetro/ (27rr|Q p^t^^) 
in units of ADUs per pixel, where fpetro and rso.Petro are the Pet- 
rosian flux and Petrosian half-light radius, respectively). For n > 4 
galaxies with a normalised sky difference less than 0.01, meaning 
that the sky difference is at most a minor issue, the n and r^o dis- 
agreements more or less reflect the systematic underestimates seen 
in the NYU-VAGC fitting of simulated galaxies. As Asky / Ipetro 
increases, the NYU-GALFIT disagreements grow and we see a 
trend of larger NYU-VAGC underestimates for galaxies with higher 
n, as expected when the local sky estimate includes more of the 
light belonging to each galaxy. 



A2 Low Sersic Galaxies 

We have outlined how the NYU-VAGC and GALFIT methods both 
do very well in fitting pure-Sersic simulations with nin < 3, there- 
fore we expect minor differences when comparing fits to real galax- 
ies with disk-like profiles. However, we find that the fit parameters 
from the two procedures differ in two ways for n < 3 galaxies, as 
shown in Figure lATI First, the NYU-VAGC fits have systematically 
higher Sersic values than the GALFIT fits, which is inconsistent 
with the results from the simulations. Second, there is a systematic 
offset of about 0.2 mag between the magnitude of NYU-VAGC fits 
and our GALFIT fits in the sense that NYU-VAGC finds fainter 
fluxes. We note that the offset appears to be independent of the dif- 
ference in the sky ( Asky = local sky — global sky) used in each 
fitting procedure, as shown by the similarity between the red (small 
Asky) and the black (large Asky) points for n < 3 galaxies in 
Figure lATI 

We suspect that the NYU-VAGC procedure of fitting a 1-D 
Sersic model to azimuthally averaged annuli overestimates Sersic 
indices for disk galaxies. We check the distribution of Sersic in- 
dices for the whole NYU-VAGC and find that the number of galax- 
ies with 0.5 < riNYU-vAGC < 1.0 is much less than those with 
1.0 < nNYU-vAGC < 1-5. Thi s results in conflic t s with other 
observations of disk galax ies (e.g. [Driver et al.ir2006l ; Ivan der Well 
l2008l :l Haussler et al 1 l2007l) . To test our suspicion, we visually in- 
spect galaxies from our sample with a late-type fit (n < 2) by 
GALFIT, but an early-type fit (riNYU-VAGC > 2.5) by NYU- 
VAGC. To exclude the sky influence, we restrict our inspection to 
36 galaxies with |Asky| < 0.5 ADU. At least two thirds of these 
galaxies have very obvious spiral features as expected for galax- 
ies with disk-dominated light profiles. Another 20 percent have 
disturbed morphologies or very bright nearby stars, which could 
cause spurious fits. A majority of the s pirals are inclin e d wit h 
b/a < 0.5. As clearly demonstrated by ISaiUn & HarrisI j2008l) . 
unyu-vagc is systematically overestimated for more inclined 
galaxies. This effect is the result of edge-on or inclined galaxies 
having steeper azimuthally averaged radial profiles because the av- 
eraged flux from the narrow outer part of such galaxy is decreased 
by being smoothed over a large circular area. 
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Figure A2. The NYU-GALFIT discrepancies shown in Figure lAll as a func- 
tion of sky difference, divided into GALFIT Sersic index bins as shown by 
the colour coding. Here the sky difference is expressed by the ratio between 
Asky = local sky — global sky and the average Petrosian surface bright- 
ness (Petrosian quantities ai'e directly drawn from SDSS tables). 

A3 Comparing Sersic Magnitude Estimates 

Besides the tendency to overestimate the Sersic indices of ac- 
tual disk-dominated galaxies when using ID fits to azimuthally- 
averaged radial profiles, we also explore the offset between the 
NYU-VAGC and the GALFIT Sersic magnitudes in more detail. 
For this exercise we use an independent measure of galaxy flux, 
the SDSS Petrosian magnitude, to anchor our comparisons of dif- 
ferent Sersic magnitudes from the two methods. The SDSS photo- 
metric pipeline calculates the flux within a circular aperture equal 
to two times the Petrosian radius, which provides an approximate 
total galaxy magnitude. It is well known, however, that Petrosian 
magnitudes systematically miss some flux whe n applied to d iffer- 
ent Sersic model profiles. As shown in Grahar net alj ( |2005|) . the 
Petrosian magnitude misses very little flux for an n = 1 profile, 
but for an n = 4 galaxy it will und erestimate the brightn ess by 
about 0.2 magnitudes. We note that the lGraham et al. calcu- 
lations are valid only when the sky is known perfectly. Any un- 
der/overestimation of the sky background will increase/decrease 
the discrepancy between th e Petrosian and the Sersic magni- 
tudes. Using the formalism o f lGrahametij] ( l2005l) . we predict the 
Petrosian-Sersic magnitude offset (Amag = Petrosian — Sersic) 
under the influence of different amounts of overestimation of the 
real sky by subtracting a range of background pedestals to each 
Sersic model before measuring the Petrosian flux. Our Amag pre- 
dictions for different n are shown in Figure I A3 1 (top panel) as a 
function of the sky overestimate (Asky), normalised by the Pet- 
rosian surface brightness (as in Figure [A2t . When the normalised 
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Figure A3. The difference between Petrosian and Sersic model quantities 
as a function of sky offset normalised to tlie Petrosian surface briglitness. 
Cases with different Sersic index are colour coded as indicated in the top 
panel. Solid lines are for our predictions based on a pure Sersic model (see 
text for details), for which Asky means the imaginary overestimation of sky 
background. The small squares are for each galaxy in our GALFIT sample, 
for which Asky means the SDSS local-global sky difference. The filled cir- 
cles with error bars show the mean and standard deviation of our sample in 
different bins. We also compare with the mean and standard deviation (tri- 
angles with error bars) of the NYU-VAGC fitting parameters, which should 
be plotted at Asky = 0, but are shifted a little to allow them to be plotted 
on a log-scale plot. 



same local sky as the Petrosian quantities, all the galaxies with 
NYU-VAGC fits have Asky = by definition. Therefore, the val- 
ues of Amag and rso.pctro/fso.Sersic for the NYU-VAGC Sersic 
results for our sample (triangles in Figure |A3l l should satisfy the 
predictions of iGraham et alj 120051) . Yet, we see that the NYU- 
VAGC Sersic results are actually underestimates, on average, com- 
pared to the predictions. For example, the n = 1 galaxies have 
Sersic magnitudes that are 0.1 m^Lg fainter than the Petrosian mea- 
surement, which is inconsistent with either the predictions or the 
definition of the two magnitudes. By definition, Sersic magnitudes 
are based on a model flux integrated to infinity, thus there is no 
reason for such a magnitude to be fainter than the Petrosian aper- 
ture magnitude, which only includes light our to some radius. It 
is unclear why the NYU-VAGC Sersic fitting procedure produces 
fainter magnitudes than expected, but this effect combined with the 
nonzero Asky values explain the systematic 0.2 mag offset that we 
find between the GALFIT and the NYU-VAGC Sersic magnitudes 
for n < 3 galaxies (Figure IaTT i. Based on the above analysis, we 
conclude that if the two assumptions of our fitting are valid, i.e. 
assuming a Sersic model is the correct model and that the SDSS 
global sky is an accurate measure of the true sky background, then 
our GALFIT fitting of the SDSS data returns more accurate mea- 
surements for the structural parameters of galaxies than those in the 
NYU-VAGC. 

This paper has been typeset from a TgX/ LSI^ file prepared by the 
author. 



sky difference is less than 10~ ^, our predictions co nverge as ex- 
pected to the values claimed in lGraham et al.l ( |2005|) . However, as 
the overestimates of the local sky increase, underestimates of the 
Petrosian magnitude for different n grow systematically. Likewise, 
we also make predictions for the offsets between the Petrosian and 
the Sersic half-light radii (r5o,pctro/''50,Scrsic) as a function of sky 
offsets and plot these in the bottom panel of Figure |A3l 

In Figure I A3 1 we also compare our GALFIT results for ac- 
tual galaxies (small squares) to the sky dependent predictions. Our 
working assumptions are: (1) a Sersic model is a reasonable model 
to describe galaxy light profiles, and (2) the SDSS global sky is 
a good measurement of the real sky and is preferred to using the 
SDSS local value. Here Asky is the local-global sky difference in 
SDSS and all the Petrosian results are measured using the local sky. 
In the upper panel, we see that our fit results are close to the Amag 
predictions for a wide range of sky differences and Sersic indices. 
We also find fair agreement between our rso results and the pre- 
dictions in the lower panel of Figure lA3l suggesting that our fitting 
results are self-consistent under the two assumptions above. 

Finally, we compare the NYU-VAGC results for the real data 
with the predictions. Given that the NYU-VAGC fitting uses the 
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